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Abstract 



Lattice QCD simulations of the potential between two baryons, each con- 
taining a heavy quark and two light quarks, such as the AqAq potential, will 
provide insight into the nucleon-nucleon interaction. As one-pion exchange 
does not contribute to the AqAq potential, the long-distance behavior is 
dominated by physics that contributes to the intermediate-range attraction 
between two nucleons. We compute the leading long-distance contributions to 
the AqAq potential in QCD and in partially-quenched QCD in the low-energy 
effective field theory. 
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I. INTRODUCTION 



The nucleon-nucleon (NN) potential has been studied for many decades, and modern 
NN potentials have been constructed to reproduce all NN scattering data below ~ 500 MeV 
to very-high precision [1]. It is well established that the long-distance component of the po- 
tential is due to one-pion exchange (OPE). A clear demonstration of this point is provided 
in work by the Nijmegen group where the charged and neutral pion masses are free param- 
eters in the OPE component that are fit to the NN phase-shift data [2]. The fit values of 
the masses reproduce with high accuracy the physical values of the pion masses. Somewhat 
more controversial are the intermediate-range and short-range components of the NN poten- 
tial. In the simplest potentials based on one-meson exchange (OME) the intermediate-range 
attraction is supplied by the exchange of the (fictitious) "a" -meson, while the short-range 
repulsion is supplied by the exchange of the p and uj mesons. In modern potentials, such as 
the AVis potential, the intermediate- and short-range components of the NN potential are 
not assumed to arise from the exchange of a single meson and ad hoc smooth functions are 
used to best reproduce the NN phase-shift data [1]. The last decade has seen considerable 
progress towards formulating NN interactions and multi-nucleon interactions in an effective 
field theory (EFT) consistent with the approximate chiral symmetry of QCD [3]. In fact, re- 
cently it was shown that an NN potential computed with the EFT better reproduces the NN 
scattering data than the potentials based on OME [4] . More precisely, the two-pion exchange 
(TPE) component of the potential provides a better description of the intermediate-range 
component than the exchange of a "a" -meson [4]. Further it was shown that the extracted 
values of coefficients of one-nucleon local operators in the chiral Lagrangian that gives rise 
to the potential are consistent with their determination in the nN sector [4]. This strongly 
suggests that the intermediate-distance components of the NN potential can be described by 
the exchange of uncorrelated pions with interactions constrained by chiral symmetry. This 
is a very encouraging situation, as ultimately one wants to make a direct and quantitative 
connection between QCD and nuclear physics using the language of effective field theory [3] . 

It is through lattice QCD simulations that definitive predictions in nuclear physics will be 
made directly from QCD. Given the current state of technology, simulations of multi-nucleon 
systems are intractable, but realistic simulations of two-nucleon systems are presently fea- 
sible. Arguably, the most promising method is to calculate scattering phase shifts directly 
using Luscher's finite-volume algorithm, which, for instance, expresses the ground-state en- 
ergy of a two-particle state as a perturbative expansion in the scattering length divided by 
the size of the box [5]. This method has been used successfully to study nn scattering [6]. 
There has also been an attempt to compute NN scattering parameters in lattice QCD; 
Ref. [7] computes the S-wave scattering lengths in quenched QCD (QQCD) using Luscher's 
method. 

A second approach to the NN system on the lattice is to study the simplified problem 
of two interacting heavy-light particles [8-13]. It has been suggested that lattice QCD 
simulations of the potential between hadrons containing a heavy quark will provide insight 
into the nature of the intermediate-range force between two nucleons [14]. In the heavy- 
quark limit, the kinetic energy of the heavy hadrons is absent and the lowest-lying energy 
eigenvalues, which can be measured on the lattice, are given by the interaction potential. 
To our knowledge, all discussions to date have addressed two heavy mesons, say BB. This 
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case is somewhat complicated by the fact that B and B* are degenerate in the heavy-quark 
limit and therefore require a coupled-channel analysis. In this work we consider the AqAq 
interaction (where the Aq is the lowest-lying baryon containing one heavy quark, Q), which 
does not suffer from this complication, and we compute the leading contribution to the 
AqAq potential in the chiral expansion. Since the Aq is an isosinglet, there is no OPE, and 
the leading large-distance behavior is governed by TPE, which is physics analogous to the 
intermediate-range attraction in the NN potential. 

Presently, unquenched lattice simulations with the physical values of the light-quark 
masses are prohibitively time-consuming, even on the fastest machines. While some 
quenched calculations can be performed with the physical quark masses there is no limit in 
which they reproduce QCD, and consequently they should be considered to be a warm-up 
exercise for the "real thing". Relatively recently it was realized that partially-quenched 
(PQ) simulations, in which the sea quarks are more massive than the valence quarks, pro- 
vide a rigorous way to determine QCD observables and are less time-consuming than their 
QCD counterparts. It is with PQ simulations that nuclear-physics observables will first be 
calculated rigorously from QCD. As PQQCD reproduces QCD only in the limit in which the 
sea-quark masses become equal to the valence-quark masses which, in turn, are set equal 
to the physical values of the light quark masses (we call this the QCD limit), there are 
some interesting features of the PQ theory that are distinct from nature away from the 
QCD limit. In QCD, the long-distance component of the NN potential is due to OPE, as 
discussed above. However, in PQQCD there is also a contribution from the exchange of 
the r^-meson (in the theory with two flavors of light quarks). In QCD such an exchange is 
suppressed due to the large mass of the rj compared to the it. However, in PQQCD the r\ 
propagator has a double-pole component that depends on the pion mass due to the hairpin 
interactions with a coefficient that depends upon the difference between the masses of the 
sea and valence quarks. Therefore, away from the QCD limit the long-distance component 
of the NN potential is dominated by one-eta exchange (OEE) and falls exponentially with 
a range determined by the pion mass, ~ e~ m7rr , as opposed to the familiar Yukawa type 
behavior [15]. This exponential behavior is even more pronounced in quenched calculations. 
Therefore, it is likely that quenched and PQ calculations of the NN potential will unfortu- 
nately be dominated by quenching and partial-quenching artifacts. While this seems like a 
negative state-of-affairs, there is hope that the PQ simulations will tell us something about 
the NN interaction. An EFT with consistent power-counting describing the multi-nucleon 
sector [16] has recently been developed, and its partially-quenched counterpart was straight- 
forwardly constructed [17]. The NN scattering amplitudes in the S-waves were computed 
out to next-to-leading order (NLO) in the EFT power-counting, while the contributions to 
the higher partial waves, resulting from OPE and OEE, were computed at leading order. 
The hope is that finite-volume lattice calculations [5] can be performed to extract the NN 
phase-shifts directly from PQQCD through the use of the scattering amplitudes computed 
in Ref. [17]. Unfortunately, given the violence that quenching and partial-quenching do to 
the long-distance component of the NN potential, there is little hope that such calculations 
can shed light on the nature of the intermediate-range component of the NN force. 
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II. THE AqAq POTENTIAL IN QCD 



The lowest-lying baryons containing a single heavy quark can be classified by the spin of 
their light degrees of freedom (dof), sj, in the heavy-quark limit, mq — > oo. Working with 
two light flavors, u and d quarks, the light dof of the isosinglet baryon, Aq, have si = 0, while 
the light dof of the isotriplet baryons, Sq 1 ' and Sq 1,0 *, (the superscript denotes the third 
component of isospin) have si = 1 1 . In the heavy-quark limit the spin-| Eg 1 ' baryons are 
degenerate with the spin-| Sq 1 ' * baryons, but are split in mass from the Aq by an amount 
that does not vanish in the chiral limit. The Lagrange density describing the interactions of 
these heavy baryons with the pions is [18,19] 

C = iK Q v ■ DA Q - iS^v ■ DS§ + Alf^Sgj + ig 2 e, l/aX T' ij v» S} k 

+ 93 (Aq^' (A„)f S% + h.c) , (1) 

where is the superfield containing both the Sq 1 ' and Sq 1 ' *, — | (£,9^ — ^d^ 
is the axial-vector field of pions, is the heavy-hadron four-velocity, (72,3 are coupling 
constants that must be determined either from data or from lattice QCD, and A is the 
Eq^ — Aq mass splitting in the heavy-quark limit. Experimentally, the measured width of 
the E* ++ is r(S* ++ ) = 17.9±|;1 ± 4.0 MeV [20], which fixes the axial coupling g 3 to be 
\gs\ = 0.95 ± 0.08 ± 0.08 in the charmed sector. Recently, it has been conjectured that in 
QCD, <73 = 1 and gi = 2 in the chiral limit [21]. The coupling g 2 might be measured through 
the radiative decays of the [22]. 

As the light dof in the Aq have s\ = in the heavy-quark limit, the light-quark axial 
current matrix element vanishes, and thus there is no AqAq7t interaction at leading order 
in the heavy quark expansion. This means that there is no OPE (or OEE) contribution to 
the AqAq potential in QCD and PQQCD, and therefore there is no long-distance compo- 
nent in the AqAq potential. It is the two-pion exchange box and crossed-box diagrams, as 
shown in Fig. 1, that provide the longest-distance interaction between two Aq's. In addi- 
tion, there are local four-AQ operators at the same order in the chiral expansion but such 
local interactions give coordinate-space delta-functions. Not only are such local interactions 
allowed, they are required to absorb the divergences arising in the box and crossed-box di- 
agrams in Fig. 1. Analogous diagrams to those in Fig. 1 in the two-nucleon sector provide 
part of the intermediate-range component of the NN potential. However, it is important 
to realize that there are additional interactions that contribute to the intermediate-range 
component of the NN potential in the chiral expansion, for instance contributions from 
the Weinberg- Tomazawa term and also from higher-dimension NNirn vertices. Therefore, 
while the AqAq potential provides a window into the nature of the intermediate-range NN 
interaction, it certainly does not provide a complete description. 

If the Aq and Sq- 1 were degenerate, we would be required to solve the coupled-channel 

system with AqAq and Eq^Eq^ coupled to / = 0. In the charmed sector the E c — A c mass 
splitting is A = 167.1 MeV and the E* - A c mass splitting is A = 232.7 MeV, and we use 



For three light flavors the baryons fall into a 6 © 3 of SU(3). 
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FIG. 1. The box and crossed-box diagrams that give the longest-distance component of the AqAq 
potential. The vertices are taken from the Lagrangian in eq. (1). Heavy-quark symmetry forbids 
AqAq intermediate states in the box and crossed-box diagrams. 

the spin- weighted average of these splittings to estimate A ~ 211 MeV in the heavy-quark 
limit. There is no symmetry reason for this mass-splitting to vanish in the chiral limit, and 
hence there is no infrared divergence that requires a coupled-channel analysis. In the power- 
counting we treat A ~ and take the Ma q , Af Sq — > oo limit in evaluating the diagrams 
in Fig. 1. With this power-counting, one can directly use the Feynman rules of Heavy- 
Baryon Chiral Perturbation Theory (HB^PT) to describe the low-momentum dynamics of 
the nucleon and A-resonance, developed by Jenkins and Manohar [23], without the need to 
resum the baryon kinetic energy term as is the case for the box and crossed-box diagrams 
in the nucleon sector. Evaluating the diagrams in Fig. 1 and then Fourier transforming 
them to position space is straightforward. Using the integral representations of Ref. [24], 
the potential between two Aq's in QCD at leading order in the chiral expansion and in the 
isospin limit is 



\/ QCD (r) 



9l 



87r 3J4 r 2 



f 

JO 



d\ 7) — X(A, to^, r) 

(A 2 + A 2 ) 2 



(2) 



where for notational simplicity we have defined the function X(A, to, r) to be 

2 



X(A,m,r) = e - 2rVX " +m ' 2 



2 A + to + 



3VA 2 + m 2 | 3 



+ (A 2 + m 2 )' 



(3) 



and /tt ~ 132 MeV. A useful feature of this integral representation is the separation between 
the A and m n dependence. The potential is shown in Fig. 2. For asymptotically large 
distances, the potential is well-approximated by 



U QCD (r) 



3 g 3 m% 



9/2 



16 ttS/2 r 5/2 A Z J4 



_ g -2m„r _|_ 



(4) 



which exhibits the expected fall off with a length scale set by twice the mass of the pion and 
where the dots represent subleading contributions in the large-distance expansion. 
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FIG. 2. The logarithm of the AqAq potential in MeV as function of the distance between the 
Aq in fm. 

III. THE AqAq POTENTIAL IN PQQCD 

The extension of the heavy-baryon sector from QCD to PQQCD is straightforward 2 . 

We use the notation and methods of Ref. [26]. To establish the representations of SU(A\2) 

(*) 

containing the Aq and Eg, we consider their interpolating fields. The interpolating field 
that contains the Aq is 



-7 
"ij 



Q 



7,c 



QTQf + Q? a Qr e abc (c 75 ) a 



a,as-\f3,b 



(5) 



while the interpolating field that contains the Eq's is 



1 1 A 4 



Q 



7,c 



QTQf - Q? a Q?' b } e abc (C 7 % 



(6) 



where Q"' a is a light-quark field operator with graded flavor index i — 1, 2, 6 where indices 
1, 2, 3, 4 are fermionic and 5, 6 are bosonic. The other indices are the color index a and Dirac 
index a. The heavy-quark field operator is Q 7,c that does not carry a graded flavor index. 
Under interchange of the light-quark flavor indices we find that 



T7 _ I \ViVj T~1 C7:^ _ / \l+ViVj C7.M 



(7) 



where rji = +1 for fermionic indices and rji = for bosonic indices. It is then straightforward 
to show, following the techniques of Ref. [26], that the Aq is embedded in a 17 dimensional 
representation of SU(4\2), %j, while the Eg" 1 are each embedded in a 19 dimensional rep- 
resentation of SU{4\2), S i:j . Under SU{2) q SU{2) j <g> SU(2)~ q the ground floor of level A 



2 The extension of the heavy-baryon sector to quenched QCD was carried out in Ref. [25] 
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of %j transforms as a (1,1,1) and contains the Aq, the first floor of level A contains four 
baryons transforming as a (2, 1, 2), while the ground floor of level B contains four baryons 
that transform as (2,2,1). Similarly, the ground floor of level A of Sij transforms as a 

(3, 1, 1) and contains the Eq^'s, the first floor of level A contains four baryons transforming 
as a (2,1,2), while the ground floor of level B contains four baryons that transform as 
(2,2,1). The remaining members of each irreducible representation do not contribute to 
the potential at the order we are working. The Lagrange density describing the dynamics 
of the lowest-lying baryons and the pseudo-Goldstone bosons is 

C = i(Tv DT) - i (Tv ■ DS^j + A 

+ i 92 Wa (S 11 v» A* S x ) + y/2 g 3 (7 A" <S M ) , (8) 
where the covariant derivatives are 

d,^ = + (y^%. + {-r^ ) {v,) ] m % m , (9) 

and similarly for the S^. The index contractions are 

Aq Z Q , Z Q A q 
« » 

M , X \ / M , X 




FIG. 3. Contributions to the crossed-box diagram from mesons and baryons containing ghost 
and sea quarks. The vertices are taken from the Lagrangian in eq. (8). The Zq are baryons 
containing one sea quark that transform as a (2,2,1) under SU{2) q ® SU(2)j <g) SU(2)g , while 
the Zq are baryons containing one ghost quark that transform as a (2, 1, 2). The x we fermionic 
mesons containing one ghost quark while the M are mesons containing one sea quark. 



(ttt)= T a > ji rt r A0 - , (3" r s,) = s a ^ s^ i3 
(tyyt)= T a ^ Yl r M , pry s„) = r i* y\ s M] 
(tty^s,) =T a ' 3l vi (y)J , (10) 

where T acts on the Dirac indices and Y — > UYW acts in flavor space. 

In PQQCD the box diagram in Fig. 1 is unmodified due to the fact that one cannot 
exchange two mesons containing either a ghost quark or a sea quark and have only two 
baryons (containing ghost and sea quarks) in the intermediate state. However, the crossed- 
box diagram is modified by quenching and partial quenching, as shown in Fig. 3. Writing 
the potential in PQQCD as V PQ (r) = V QCD (r) + SV(r), we find that 

o 4 f°° A 2 - A 2 
WW = +M$Jy J, d\ (l(X,m S v,r) - J(A, m,, r) ) , (11) 
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FIG. 4. The left panel shows I/ p Q(r) evaluated at r = 1 fm as a function of the meson mass 
m SV; while the right panels shows y P( ^(r) evaluated at r = 2 fm. The vertical axis is in units of 
MeV. When msv = m^ the value ofV P( ^(r) is equal to V® CD (r). 



where msv is the mass of a meson containing one sea quark and one valence quark. In Fig. 4 
we show I/ p Q(r) evaluated at r = 1 fm and r = 2 fm as a function of msv- I n contrast 
to the NN potential that is dominated at long-distances by partial-quenching effects, the 
AqAq potential receives comparable contributions from both the valence quarks and the 
ghost quarks, and in fact the (partial) quenching contribution is somewhat smaller than the 
QCD contribution. When msv — m n we recover the value of the potential in QCD, as has 
to be the case. 




FIG. 5. The leading hairpin contribution to the AqAq potential. The crosses denote double-pole 
rj propagators and the boxes are vertices from Lagrange density in eq. (12) 



Given the large hairpin contribution to the NN potential, it is important to estimate 
the size of the hairpin contribution to the AqAq potential. The Lagrange density for the 
leading 7] (hairpin) interactions, neglecting operators involving insertions of the light-quark 
mass matrix, is 



1 

x 



Ci (Vr/) 2 + c 2 (<9 ?7) 2 



AnAc 



(12) 



where A x is the chiral-symmetry breaking scale and are unknown constants whose natural 
size is C12 ~ 1 by naive dimensional analysis. One finds that two-insertions of the interactions 
in eq. 12, as shown in Fig. 5, give a contribution to the AqAq potential of 



S 



V 



HP, 



ss 



32tt 3 A6 



™£f r dX e ~2r^t r f 2 + A 2 + m 2\ c 2 + 2A 2 ^ 

A° jo L \r z / 



A' 1 



A 2 + 



(13) 



The remaining integral can be performed analytically in terms of modified Bessel functions 
and Struve functions, but it is less than illuminating and we do not present it here. This 
contribution is suppressed compared to the leading non-hairpin contributions as dictated 
by the EFT power counting. Numerically, we find that the hairpin contribution to the 
potential at r = 1 fm is V RP (1 fm) ~ -1 keV, while at r = 2 fm is V up (2 fm) ~ -0.06 keV 
for ci = c 2 = 1, A x ~ 1 GeV and m ss = 600 MeV. At asymptotically large distances the 
potential becomes 



V (r) 



[m 



ss 



m 2 ) 2 ml/ 2 



-2m 7T r 



64 7T 5 /2 A<3 



+ 



(14) 



where the dots represent subleading contributions in the large-distance expansion. While 
at asymptotically large distances this contribution is larger than that from the box and 
crossed-box diagrams, asymptopia finally sets in at distances at which all contributions are 
numerically insignificant. 



IV. CONCLUSIONS 



In this work we have computed the potential between two Aq's at leading order in 
effective field theory in both QCD and PQQCD. Further, we have estimated the size of 
the leading contribution from hairpin interactions in PQQCD. We find that the partially- 
quenched AqAq potential does not suffer from some of the (partial-) quenching problems 
that plague the NN potential due to the absence of single pseudo-Goldstone exchange. The 
potentials we have computed in this work will allow for the chiral extrapolation of lattice 
calculations performed with unphysically large sea quark masses. 

As these potentials fall off with a mass scale set by ~ 2m n , they are quite small for baryon 
separations greater than r ~ 1.5 fm. Therefore, the theoretical advantages of studying this 
system to learn about the NN potential may be undermined by the difficulties in extracting 
a signal from lattice simulations. However, the simplifications introduced by only having two 
light quarks, and a single infinitely-massive quark to fix the inter-baryon separation makes 
this system a prime candidate for studying inter-baryon interactions. We look forward to 
the Jlab/MIT group [14] carrying out their program in this area. 
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